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Abstract 

Using a new representation for the quasiclassical Green function of the Dirac 
equation in a Coulomb held, analytical expressions for the high-energy small-angle 
Delbriick scattering amplitudes are obtained exactly in the parameter Za. Magni- 
tudes of the amplitudes coincide with the previous results. However, the structure of 
the expressions obtained is much more simple, which considerably facilitates numer- 
ical calculations. 



The coherent scattering of photons in the electric fields of atoms via virtual electron- 
positron pairs (Delbriick scattering |1[) has been intensively investigated both theoretically 
and experimentally . From the theoretical point of view this process is interesting since 
the higher orders of the perturbation theory in the parameter Za (Z\e\ is the nucleus 
charge, a = e 2 = 1/137 is the fine-structure constant, e is the electron charge, h = c = 1) 
drastically change the amplitudes as compared to the first (Born) approximation. Delbriick 
scattering at high photon energy uo ^> m (m is the electron mass) is the most favorable to 
compare the theoretical and experimental results. Recently, the most accurate Delbriick 
scattering cross section measurements have been performed in the Budker Institute of 
Nuclear Physics || at photon energies 140 — 450 MeV and scattering angles 2.6 — 16.6 
mrad. The Delbriick scattering amplitudes exact in the parameter Za at high photon 
energies and small scattering angle were obtained in j|, |5|, || by summing in a definite 
approximation the Feynman diagrams with an arbitrary number of photons exchanged 
with a Coulomb center. Another representation of these amplitudes was derived in |7], [8] 
using the quasiclassical Green function of the Dirac equation in a Coulomb field. Such a 
Green function in an arbitrary spherically symmetrical external field was obtained in |], |I~0 
where Delbriick scattering in a screened Coulomb field was investigated as well. 



In recent papers [12], [13], |14| the new approach was developed, which significantly sim- 
plifies the consideration of high-energy processes in a spherically symmetrical external field. 
It was used to obtain the exact in Za photon splitting amplitudes. These results are used 
in the data processing of the first successful experiment on the high-energy photon split- 
ting in the electric fields of atoms recently performed in the Budker INP (the preliminary 
experimental results are published in [jTTH) . 



In the present paper the approach developed in [12|, [13], [I4J] is applied to the calculation 
of the high-energy small-angle Delbriick scattering amplitudes. Numerically, the new am- 
plitudes coincide with the previously known results. However, the numerical calculations 
by means of the new formulas, having actually only double integral, are significantly easier. 

As shown in [[|, it is convenient to express the Delbriick scattering amplitude via the 
Green function of the squared Dirac equation D(r, r'\ e): 

D(r 1 ,r 2 | £ ) = (r 1 |l/(P 2 -m 2 + zO)|r 2 ) , 

where V = 7°(e + Za/r) - 7 p, p = — iV . 

In terms of the function D(r, r'\e) the amplitude of the process reads (Pj): 
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where e±, k\ (e 2 , k 2 ) are the polarization vector and 4-momentum of the initial (final) 
photon, p 12 = — iVi j2 . By definition, the Delbriick scattering amplitude ([[]) should vanish 
at Z = 0. Therefore, one should subtract from the integrand in (fj) its value at Z — 0. We 



perform this subtraction in the explicit form below. At u ^> m the main contribution to 
the cross section comes from the momentum transfer region A ~ m corresponding to the 
small scattering angles. Then we can neglect the last term in (H) since it depends only on 
the momentum transfer A = k 2 — ki, while the amplitude at uj ^> A is proportional to uj 

(see, e.g., @). 

According to the uncertainty relation, the lifetime of the virtual electron-positron pair is 
t ~ |r2 — ri| ~ uj/(m 2 + A 2 ) and p ~ 1/A for the characteristic impact parameter. For uj ^> 
A 3> m 2 /uj the main contribution to the integral in ([I]) is given by small angles between the 
vectors r 2 , — r 1; and ki. Then the characteristic angular momentum is I ~ ujp ~ uj/ A 3> 1 
and the quasiclassical approximation is valid. Note that the screening should be taken 



into account only at A 
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<C m, where r c is the screening radius 
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in the Thomas- Fermi model). In the present paper we consider the momentum transfer 
region A 3> m 2 /uj, r" 1 , which gives the main contribution to the total cross section of the 
process. We emphasize that the contribution of higher orders of perturbation theory in Za 
(Coulomb corrections) is not affected by screening and the corresponding expressions are 
valid for any A <C uj. The modification of the Born contribution at A ~ m 2 juj was studied 
in , the effects of screening were discussed in detail in . 

In papers |T2], [Tj|, the convenient representation for the quasiclassical Green function 
-D(i"i, r 2 |e) of the squared Dirac equation in a Coulomb field was obtained. For small angles 
between r 2 , —X\ and the z axis, we have 
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where a = 7°7 , k 2 = e 2 — m 2 , A = e/k , q, 6\, 6 2 are two-dimensional vectors in the xy 
plane, Q\ = ri±/r x , 6 2 = r 2 ±/r 2 . Expression (Q) contains only elementary functions, and 
the angles #i 6 2 appear only in the factor exp[iq(6i + 2 )]. Therefore, representation (|2|) 
for the Green function is very convenient for calculations. 

We direct the z axis along ki. The z component of the polarization vector e 2 can be 
eliminated owing to the relation e 2 k 2 = 0, which leads to (e 2 ) z = —e 2 A./uj . After that, 
within the small-angle approximation one can neglect the difference between the vector 
(e 2 )j_ and the polarization vector of a photon, propagating along the z-axis and having the 
same helicity. So, the amplitudes of the process are expressed via the transverse polarization 
vectors e and e*, corresponding to the positive and negative helicities, respectively. It is 
sufficient to calculate two amplitudes, for instance, M ++ and M + _. Other amplitudes 
(M and M_ + ) can be obtained by the substitution e ^ e*. 

Substituting (fj) into (|), we expand the amplitudes at small angles when dridr 2 ps 
r 2 r 2 dridr 2 d0 1 d0 2 . Taking the trace and performing the elementary integration over the 
angles Ox and B 2 , we obtain 
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Here 
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the function T for different polarizations has the form 
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and following notation is introduced k = uj — e, Q = qi + q 2 q = qi — q 2 . When deriving 
(0) , we integrate by parts over qi, q 2 so that the integrand contains the parameter Za 



only in the factor [(qi/q$ 



\2iZc 



1]. Besides, we make the substitution r 12 — > (sk/uj)^?. 



The expression T ++ in (|5p can be simplified if we pass temporarily in (^) from the variables 
ri, r 2 to i? = T\T%j{r\ + r 2 ) , £ = ri/r 2 and integrate by parts over i? the term proportional 
to Q 2 in square brackets in (|) . As a result, we get 
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Let us pass from the variables qi, q 2 to q and Q. Then the integral over q acquires the 
form: 
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As shown in [Q, this integral is equal to 
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Using this representation and the parametrization 
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where x is a two-dimensional vector lying in the same plane as Q, we take the integrals in 
(H) first over r 1; then over Q and over r 2 . Finally, we have 
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where g± = |q±|, q± = q±A, v = x + q/2 + A [e — n)/2u. 
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For the further integration it is convenient to rewrite the expression (10) in another 
form. Using the identities 
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and integrating by parts over q, we find 



M 4 
M 4 



2a(Za)m' 

TT 2 A 2 UJ 2 



u / \ 2iZa 

J* Hz) /■ 



dx. 



x 2 + m 2 ) 2 (v 2 + m 2 ) 



x 



X 



q+ _ q- 



u 2 x. — (e 2 + k 2 )\ 
AeK,(e\)e 



;n) 



Applying the Feynman parametrization for the denominators, we take the integral over x 
and pass from the variable e to s = 2e/uo — 1. We obtain 
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We perform the integration over q by means of a trick used in |13| , [14]]. Let us multiply the 
integrand in (0) by 

= (q 2 + A 2 ) f ^((q - A/y) 2 - A 2 (l/y 2 - 1)), 

- 7 - 1 \y\ 

change the order of integration over q and y, and make the shift q — > q + A/y. After 

that the integration over q can be done easily. By changing the variables y = tanh(r — To), 

where 
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we express the integral over r in terms of two functions, the same as in |14 
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where a 2 = AB/[(B + (1 + s) 2 )(-B + (1 — s) 2 )], i^(x) is the derivative of the Legendre 
function. Note that a 2 < 1 for any s and B > 0. 



The final result for the Delbriick scattering amplitudes reads: 
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Let us derive now the asymptotics of the amplitudes flTBl) at A < m and A ^> m. 
In the case A < m, when 5 ~ m 2 /A 2 » 1, o 2 « 4/5 < 1, and r ~ 1/B < 1 , the 
calculations are especially easy. Then the functions T\p, have the asymptotic form: 

fi«l, ^2 = -2Zaa 2 [lna + C + Re^(l + ^a)], 

where C = 0.577. . . is the Euler constant, ip(x) = dlnT(x)/dx. Substituting the asymp- 
totics of the functions T\p into (|I5D, we find 
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To obtain the asymptotics at A > m, it is convenient to start from (|T2"|). The main 
contribution to the integral over t comes from the region 1 — t ~ m 2 /A 2 <C 1. Taking in 
(O) the integral over £ in this approximation, we have: 



M 4 



a(Za)cj 
A 2- 



7T 



ds 



dq 



9+ 



2iZc 



q + sA) 2 yg_ 



;i + s 2 )(q + sA) 



q+ q_\ I 

g 2 . g 2 ii 2(l- S 2 )(e,q + S A)e 



18) 



Substituting (|H|) into this formula, we take the elementary integral over q, and then over 
s. After that we perform the substitution y = tanhr and obtain at A ^> m: 
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The asymptotics ([HP and (|T9|) coincide with the results of || §]. Additionally, we checked 
numerically for different Z that the magnitudes of the amplitudes (|15|), as should be, 
coincide with the results of || [F], || at intermediate values of the momentum transfer A. 
The expression (|T5| ) for the Delbriick scattering amplitudes is a double integral, being 
essentially simpler than the known representations || ^, [|. Therefore, the formula C [3~5|) is 
very convenient for numerical calculations. 

The calculation of the Delbriick scattering amplitudes, performed in the present paper, 
demonstrate again the advantage of our approach in the consideration of different high- 
energy QED processes in a Coulomb field. 
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